Fluid suppression in amide proton transfer-weighted (APTw) CEST imaging: New theoretical insights and clinical benefits.

PURPOSE: Amide proton transfer-weighted (APTw) MRI at 3T provides a unique contrast for brain tumor imaging. However, APTw imaging suffers from hyperintensities in liquid compartments such as cystic or necrotic structures and provides a distorted APTw signal intensity. Recently, it has been shown that heuristically motivated fluid suppression can remove such artifacts and significantly improve the readability of APTw imaging. THEORY AND METHODS: In this work, we show that the fluid suppression can actually be understood by the known concept of spillover dilution, which itself can be derived from the Bloch-McConnell equations in comparison to the heuristic approach. Therefore, we derive a novel post-processing formula that efficiently removes fluid artifact, and explains previous approaches. We demonstrate the utility of this APTw assessment in silico, in vitro, and in vivo in brain tumor patients acquired at MR scanners from different vendors. RESULTS: Our results show a reduction of the CEST signals from fluid environments while keeping the APTw-CEST signal intensity almost unchanged for semi-solid tissue structures such as the contralateral normal appearing white matter. This further allows us to use the same color bar settings as for conventional APTw imaging. CONCLUSION: Fluid suppression has considerable value in improving the readability of APTw maps in the neuro-oncological field. In this work, we derive a novel post-processing formula from the underlying Bloch-McConnell equations that efficiently removes fluid artifact, and explains previous approaches which justify the derivation of this metric from a theoretical point of view, to reassure the scientific and medical field about its use.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

PURPOSE: To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction. METHODS: Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF), R 2 * $$ {\mathrm{R}}_2^{\ast } $$ , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and between R 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM in in vivo study. RESULTS: Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.31 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.05, R 2 $$ {R}^2 $$ = 0.999) and Cartesian without motion ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.32 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.04, R 2 $$ {R}^2 $$ = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.00261 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.524, R 2 $$ {R}^2 $$ = 0.977) compared to motion-averaged reconstruction ( QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.0021 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.572, R 2 $$ {R}^2 $$ = 0.723) in ROI-based linear regression analysis between R 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM. CONCLUSION: Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

GRASPNET: Fast spatiotemporal deep learning reconstruction of golden-angle radial data for free-breathing dynamic contrast-enhanced magnetic resonance imaging.

The purpose of the current study was to develop a deep learning technique called Golden-angle RAdial Sparse Parallel Network (GRASPnet) for fast reconstruction of dynamic contrast-enhanced 4D MRI acquired with golden-angle radial k-space trajectories. GRASPnet operates in the image-time space and does not use explicit data consistency to minimize the reconstruction time. Three different network architectures were developed: (1) GRASPnet-2D: 2D convolutional kernels (x,y) and coil and contrast dimensions collapsed into a single combined dimension; (2) GRASPnet-3D: 3D kernels (x,y,t); and (3) GRASPnet-2D + time: two 3D kernels to first exploit spatial correlations (x,y,1) followed by temporal correlations (1,1,t). The networks were trained using iterative GRASP reconstruction as the reference. Free-breathing 3D abdominal imaging with contrast injection was performed on 33 patients with liver lesions using a T1-weighted golden-angle stack-of-stars pulse sequence. Ten datasets were used for testing. The three GRASPnet architectures were compared with iterative GRASP results using quantitative and qualitative analysis, including impressions from two body radiologists. The three GRASPnet techniques reduced the reconstruction time to about 13 s with similar results with respect to iterative GRASP. Among the GRASPnet techniques, GRASPnet-2D + time compared favorably in the quantitative analysis. Spatiotemporal deep learning enables reconstruction of dynamic 4D contrast-enhanced images in a few seconds, which would facilitate translation to clinical practice of compressed sensing methods that are currently limited by long reconstruction times.

Dynamic contrast-enhanced MRI parametric mapping using high spatiotemporal resolution Golden-angle RAdial Sparse Parallel MRI and iterative joint estimation of the arterial input function and pharmacokinetic parameters.

The aim of this work is to develop a data-driven quantitative dynamic contrast-enhanced (DCE) MRI technique using Golden-angle RAdial Sparse Parallel (GRASP) MRI with high spatial resolution and high flexible temporal resolution and pharmacokinetic (PK) analysis with an arterial input function (AIF) estimated directly from the data obtained from each patient. DCE-MRI was performed on 13 patients with gynecological malignancy using a 3-T MRI scanner with a single continuous golden-angle stack-of-stars acquisition and image reconstruction with two temporal resolutions, by exploiting a unique feature in GRASP that reconstructs acquired data with user-defined temporal resolution. Joint estimation of the AIF (both AIF shape and delay) and PK parameters was performed with an iterative algorithm that alternates between AIF and PK estimation. Computer simulations were performed to determine the accuracy (expressed as percentage error [PE]) and precision of the estimated parameters. PK parameters (volume transfer constant [Ktrans ], fractional volume of the extravascular extracellular space [ve ], and blood plasma volume fraction [vp ]) and normalized root-mean-square error [nRMSE] (%) of the fitting errors for the tumor contrast kinetic data were measured both with population-averaged and data-driven AIFs. On patient data, the Wilcoxon signed-rank test was performed to compare nRMSE. Simulations demonstrated that GRASP image reconstruction with a temporal resolution of 1 s/frame for AIF estimation and 5 s/frame for PK analysis resulted in an absolute PE of less than 5% in the estimation of Ktrans and ve , and less than 11% in the estimation of vp . The nRMSE (mean ± SD) for the dual temporal resolution image reconstruction and data-driven AIF was 0.16 ± 0.04 compared with 0.27 ± 0.10 (p < 0.001) with 1 s/frame using population-averaged AIF, and 0.23 ± 0.07 with 5 s/frame using population-averaged AIF (p < 0.001). We conclude that DCE-MRI data acquired and reconstructed with the GRASP technique at dual temporal resolution can successfully be applied to jointly estimate the AIF and PK parameters from a single acquisition resulting in data-driven AIFs and voxelwise PK parametric maps.

QQ-NET - using deep learning to solve quantitative susceptibility mapping and quantitative blood oxygen level dependent magnitude (QSM+qBOLD or QQ) based oxygen extraction fraction (OEF) mapping.

PURPOSE: To improve accuracy and speed of quantitative susceptibility mapping plus quantitative blood oxygen level-dependent magnitude (QSM+qBOLD or QQ) -based oxygen extraction fraction (OEF) mapping using a deep neural network (QQ-NET). METHODS: The 3D multi-echo gradient echo images were acquired in 34 ischemic stroke patients and 4 healthy subjects. Arterial spin labeling and diffusion weighted imaging (DWI) were also performed in the patients. NET was developed to solve the QQ model inversion problem based on Unet. QQ-based OEF maps were reconstructed with previously introduced temporal clustering, tissue composition, and total variation (CCTV) and NET. The results were compared in simulation, ischemic stroke patients, and healthy subjects using a two-sample Kolmogorov-Smirnov test. RESULTS: In the simulation, QQ-NET provided more accurate and precise OEF maps than QQ-CCTV with 150 times faster reconstruction speed. In the subacute stroke patients, OEF from QQ-NET had greater contrast-to-noise ratio (CNR) between DWI-defined lesions and their unaffected contralateral normal tissue than with QQ-CCTV: 1.9 ± 1.3 vs 6.6 ± 10.7 (p = 0.03). In healthy subjects, both QQ-CCTV and QQ-NET provided uniform OEF maps. CONCLUSION: QQ-NET improves the accuracy of QQ-based OEF with faster reconstruction.

Integrated quantitative susceptibility and R2 * mapping for evaluation of liver fibrosis: An ex vivo feasibility study.

To develop a method for noninvasive evaluation of liver fibrosis, we investigated the differential sensitivities of quantitative susceptibility mapping (QSM) and R2 * mapping using corrections for the effects of liver iron. Liver fibrosis is characterized by excessive accumulation of collagen and other extracellular matrix proteins. While collagen increases R2 * relaxation, measures of R2 * for fibrosis are confounded by liver iron, which may be present in the liver over a wide range of concentrations. The diamagnetic collagen contribution to susceptibility values measured by QSM is much less than the contribution of highly paramagnetic iron. In 19 ex vivo liver explants with and without fibrosis, QSM (χ), R2 * and proton density fat fraction (PDFF) maps were constructed from multiecho gradient-recalled echo (mGRE) sequence acquisition at 3 T. Median parameter values were recorded and differences between the MRI parameters in nonfibrotic vs. advanced fibrotic/cirrhotic samples were evaluated using Mann-Whitney U tests and receiver operating characteristic analyses. Logistic regression with stepwise feature selection was employed to evaluate the utility of combined MRI measurements for detection of fibrosis. Median R2 * increased in fibrotic vs. nonfibrotic liver samples (P = .041), while differences in χ and PDFF were nonsignificant (P = .545 and P = .395, respectively). Logistic regression identified the combination of χ and R2 * significant for fibrosis detection (logit [prediction] = -8.45 + 0.23 R2 * - 28.8 χ). For this classifier, a highly significant difference between nonfibrotic vs. advanced fibrotic/cirrhotic samples was observed (P = .002). The model exhibited an AUC of 0.909 (P = .003) for detection of advanced fibrosis/cirrhosis, which was substantially higher compared with AUCs of the individual parameters (AUC 0.591-0.784). An integrated QSM and R2 * analysis of mGRE 3 T imaging data is promising for noninvasive diagnostic assessment of liver fibrosis.

Deep neural network for water/fat separation: Supervised training, unsupervised training, and no training.

PURPOSE: To use a deep neural network (DNN) for solving the optimization problem of water/fat separation and to compare supervised and unsupervised training. METHODS: The current T2∗ -IDEAL algorithm for solving water/fat separation is dependent on initialization. Recently, DNN has been proposed to solve water/fat separation without the need for suitable initialization. However, this approach requires supervised training of DNN using the reference water/fat separation images. Here we propose 2 novel DNN water/fat separation methods: 1) unsupervised training of DNN (UTD) using the physical forward problem as the cost function during training, and 2) no training of DNN using physical cost and backpropagation to directly reconstruct a single dataset. The supervised training of DNN, unsupervised training of DNN, and no training of DNN methods were compared with the reference T2∗ -IDEAL. RESULTS: All DNN methods generated consistent water/fat separation results that agreed well with T2∗ -IDEAL under proper initialization. CONCLUSION: The water/fat separation problem can be solved using unsupervised deep neural networks.

Quantitative susceptibility mapping of the spine using in-phase echoes to initialize inhomogeneous field and R2* for the nonconvex optimization problem of fat-water separation.

Quantitative susceptibility mapping (QSM) of human spinal vertebrae from a multi-echo gradient-echo (GRE) sequence is challenging, because comparable amounts of fat and water in the vertebrae make it difficult to solve the nonconvex optimization problem of fat-water separation (R2*-IDEAL) for estimating the magnetic field induced by tissue susceptibility. We present an in-phase (IP) echo initialization of R2*-IDEAL for QSM in the spinal vertebrae. Ten healthy human subjects were recruited for spine MRI. A 3D multi-echo GRE sequence was implemented to acquire out-phase and IP echoes. For the IP method, the R2* and field maps estimated by separately fitting the magnitude and phase of IP echoes were used to initialize gradient search R2*-IDEAL to obtain final R2*, field, water, and fat maps, and the final field map was used to generate QSM. The IP method was compared with the existing Zero method (initializing the field to zero), VARPRO-GC (variable projection using graphcuts but still initializing the field to zero), and SPURS (simultaneous phase unwrapping and removal of chemical shift using graphcuts for initialization) on both simulation and in vivo data. The single peak fat model was also compared with the multi-peak fat model. There was no substantial difference on QSM between the single peak and multi-peak fat models, but there were marked differences among different initialization methods. The simulations demonstrated that IP provided the lowest error in the field map. Compared to Zero, VARPRO-GC and SPURS, the proposed IP method provided substantially improved spine QSM in all 10 subjects.

Evaluation of diffusion kurtosis imaging in stratification of nonalcoholic fatty liver disease and early diagnosis of nonalcoholic steatohepatitis in a rabbit model.

PURPOSE: To examine the feasibility of MR diffusion kurtosis imaging (DKI) for characterizing nonalcoholic fatty liver disease (NAFLD) and diagnosing nonalcoholic steatohepatitis (NASH). METHODS: Thirty-two rabbits on high fat diet with different severities of NAFLD were imaged at 3 T MR including diffusion weighted imaging (DWI) and DKI using b values of 0, 400, 800 s/mm2 with 15 diffusion directions at each b value. Apparent diffusion coefficient (ADC) was derived from the linear exponential DWI model. Mean diffusion (MD) and mean kurtosis (MK) were derived from the quadratic exponential model of DKI. Correlations between MR parameters and hepatic pathology determined by the NAFLD activity scoring system were analyzed by Spearman rank correlation analysis. Receiver operating characteristic analyses were applied to determine the cutoff values of MD, MK as well as ADC in distinguishing NASH from non-NASH. The diagnostic efficacies of MD and MK in detecting NASH were compared with that of ADC. RESULTS: Values for ADC and MD significantly decreased as the severity of NAFLD increased (ρ = -0.529, -0.904, respectively; P < 0.05). MK values significantly increased as the severity of NAFLD increased (ρ = 0.761; P < 0.05). In addition, both MD and MK values were significantly different between borderline NASH and NASH groups (MD: 1.729 ±â€¯0.144 vs. 1.458 ±â€¯0.240[×10-3 mm2/s]; MK: 1.096 ±â€¯0.079 vs. 1.237 ±â€¯0.180; P < 0.05). Moreover, there was a significantly higher area under the curve (AUC) for both MD (0.955) and MK (0.905), as compared to ADC (0.736). CONCLUSION: Diffusion kurtosis imaging was feasible for stratifying NAFLD, and more accurately discriminated NASH from non-NASH when compared with DWI.

Measurement of the liver iron concentration in transfusional iron overload by MRI R2* and by high-transition-temperature superconducting magnetic susceptometry.

PURPOSE: To compare measurement of the liver iron concentration in patients with transfusional iron overload by magnetic resonance imaging (MRI), using R2*, and by magnetic susceptometry, using a new high-transitiontemperature (high-Tc; operating at 77 K, cooled by liquid nitrogen) superconducting magnetic susceptometer. METHODS: In 28 patients with transfusional iron overload, 43 measurements of the liver iron concentration were made by both R2* and high-Tc magnetic susceptometry. RESULTS: Measurements of the liver iron concentration by R2* and high-Tc magnetic susceptometry were significantly correlated when comparing all patients (Pearson's r = 0.91, p < 0.0001) and those with results by susceptometry >7 mg Fe/g liver, dry weight (r = 0.93, p = 0.006). In lower ranges of liver iron, no significant correlations between the two methods were found (0 to <3.2 mg Fe/g liver, dry weight: r = 0.2, p = 0.37; 3.2 to 7 mg Fe/g liver, dry weight: r = 0.41; p = 0.14). CONCLUSION: The lack of linear correlation between R2* and magnetic susceptibility measurements of the liver iron concentration with minimal or modest iron overload may be due to the effects of fibrosis and other cellular pathology that interfere with R2* but do not appreciably alter magnetic susceptibility.

Rapid automated liver quantitative susceptibility mapping.

BACKGROUND: Accurate measurement of the liver iron concentration (LIC) is needed to guide iron-chelating therapy for patients with transfusional iron overload. In this work, we investigate the feasibility of automated quantitative susceptibility mapping (QSM) to measure the LIC. PURPOSE: To develop a rapid, robust, and automated liver QSM for clinical practice. STUDY TYPE: Prospective. POPULATION: 13 healthy subjects and 22 patients. FIELD STRENGTH/SEQUENCES: 1.5 T and 3 T/3D multiecho gradient-recalled echo (GRE) sequence. ASSESSMENT: Data were acquired using a 3D GRE sequence with an out-of-phase echo spacing with respect to each other. All odd echoes that were in-phase (IP) were used to initialize the fat-water separation and field estimation (T2 *-IDEAL) before performing QSM. Liver QSM was generated through an automated pipeline without manual intervention. This IP echo-based initialization method was compared with an existing graph cuts initialization method (simultaneous phase unwrapping and removal of chemical shift, SPURS) in healthy subjects (n = 5). Reproducibility was assessed over four scanners at two field strengths from two manufacturers using healthy subjects (n = 8). Clinical feasibility was evaluated in patients (n = 22). STATISTICAL TESTS: IP and SPURS initialization methods in both healthy subjects and patients were compared using paired t-test and linear regression analysis to assess processing time and region of interest (ROI) measurements. Reproducibility of QSM, R2 *, and proton density fat fraction (PDFF) among the four different scanners was assessed using linear regression, Bland-Altman analysis, and the intraclass correlation coefficient (ICC). RESULTS: Liver QSM using the IP method was found to be ~5.5 times faster than SPURS (P < 0.05) in initializing T2 *-IDEAL with similar outputs. Liver QSM using the IP method were reproducibly generated in all four scanners (average coefficient of determination 0.95, average slope 0.90, average bias 0.002 ppm, 95% limits of agreement between -0.06 to 0.07 ppm, ICC 0.97). DATA CONCLUSION: Use of IP echo-based initialization enables robust water/fat separation and field estimation for automated, rapid, and reproducible liver QSM for clinical applications. LEVEL OF EVIDENCE: 1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:725-732.

Vastly accelerated linear least-squares fitting with numerical optimization for dual-input delay-compensated quantitative liver perfusion mapping.

PURPOSE: To propose an efficient algorithm to perform dual input compartment modeling for generating perfusion maps in the liver. METHODS: We implemented whole field-of-view linear least squares (LLS) to fit a delay-compensated dual-input single-compartment model to very high temporal resolution (four frames per second) contrast-enhanced 3D liver data, to calculate kinetic parameter maps. Using simulated data and experimental data in healthy subjects and patients, whole-field LLS was compared with the conventional voxel-wise nonlinear least-squares (NLLS) approach in terms of accuracy, performance, and computation time. RESULTS: Simulations showed good agreement between LLS and NLLS for a range of kinetic parameters. The whole-field LLS method allowed generating liver perfusion maps approximately 160-fold faster than voxel-wise NLLS, while obtaining similar perfusion parameters. CONCLUSIONS: Delay-compensated dual-input liver perfusion analysis using whole-field LLS allows generating perfusion maps with a considerable speedup compared with conventional voxel-wise NLLS fitting. Magn Reson Med 79:2415-2421, 2018. © 2017 International Society for Magnetic Resonance in Medicine.